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ABSTRACT 

We review the theory and observations related to the "supcrhump" precession of ec- 
centric accretion discs in close binary sytems. We agree with earlier work, although for 
different reasons, that the discrepancy between observation and dynamical theory im- 
plies that the effect of pressure in the disc cannot be neglected. We extend earlier work 
that investigates this effect to include the correct expression for the radius at which 
resonant orbits occur. Using analytic expressions for the accretion disc structure, we 
derive a relationship between the period excess and mass-ratio with the pressure effects 
included. This is compared to the observed data, recently derived results for detailed 
integration of the disc equations and the equivalent empirically derived relations and 
used to predict values for the mass ratio based on measured values of the period excess 
for 88 systems. 

Key words: stars: binaries: close - stars: novae: cataclysmic variables - stars: dwarf 
novae - accretion, accretion discs 



1 INTRODUCTION 

Cataclysmic Variables (CVs) are a class of close binary sys- 
tem, with a typical orbital period of a few hours, where 
a white dwarf primary accretes material from its compan- 
ion via Roche lobe overflow. The secondary is normally a 
late-type main sequence star although examples with an 
evolved companion do exist (eg. GK Per). In the absence 
of a significant white dwarf magnetic field, material arrives 
at the primary after processing though an accretion disc. 
The dwarf nova subgroup show outbursts where the lumi- 
nosity of the system increases by around 2-5 mag. Although 
not strictily periodic, these recur on a typical timescale for 
each system ranging from tens of days to tens of years. The 
dwarf novae are further subdivided based on the properties 
of the outbursts. We will be interested in the SU UMa type 
where the systems show occassional superoutbursts which 
have a brighter maximum (~ 0.7 mag.) and longer duration 
(~ 5 times) than normal outbursts. 

The most favoured explanation for dwarf nova outbursts 
involves an ionization instability where the accretion disc 
mass increases until a critical surface density 
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is reached at some radius r. When this occurs, the disc 
switches into a "hot" state with higher viscosity that causes 
a larger mass transport rate through the disc and increased 
luminosity. The disc mass now steadily decreases until a sec- 
ond critical surface density 



is reached at some point. Upon meeting this condition, the 
disc transfers back to the "cold" quiescent state and the cy- 
cle repeats. In these expressions ac and an are the Shakura- 
Sunyaev viscosity parameters in the cold and hot states re- 
spective ly and the primary mass M\ is me asured in solar 
masses dCannizzo. Shatter fc Wheeler lll98l) . 

SU UMa superoutbursts also have the property of show- 
ing superhumps. Here, an additional periodity (P s h) a few 
percent longer than the orbital period (P or b) is apparent in 
the lightcurve. This is believed to arise from a precessing, 
eccentric accretion disc driven by a resonance between the 
orbiting disc material and the secondary. CVs in general 
an d the SU UMa systems in particular are well-reviewed 
by IWarnerl lll995l) . Some systems other than dwarf novae 
also show superhumps which, by analogy, are believed to 
share a common origin with an eccentric disc now per- 
manently present givin g rise to "permanent superhumps" 
jRetter fc Navlorl l200dh Similarly, some Low-Mass X-ray 
Binaries (LMXBs); analogues to CVs where the primary is 
now a neutron star, hav e also been discov ered to have su- 
perhumps eg. KV UMa dZurita et alj l2002). There is also a 
related phenomenon of "negative" superhumps which occur 
on a period a few percent shorter than P or b- This is believed 
to arise from precession of a disc warp and we will not con- 
sider these in this paper. 

The intent of this paper is to test our understanding of 
the (positive) superhump phenomenon by comparing obser- 
vation to the predictions of theoretical expressions for P s h. 
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Since this is directly observable, if we can relate it to the 
funadmental parameters of a system, we will have a method 
to indirectly measure such parameters. 



2 THEORETICAL BACKGROUND 

Lubo w| (Il991allbl Il992fl derived the final precession rate u) 
for an eccentric disc as the sum of three terms: 

<J-> = k'dyn + Wprcss + ^tran (3) 

where uja yn is the dynamical precession frequency, u prcss is 
a pressure related term and tJtran is a transient term. This 
latter contribution is related to the time derivative of the 
mode giving rise to the dynamical precession. Thus, it may 
be important in the development phase of the superhumps 
but not in steady state. As a result, we shall not consider it 
in detail except to note that it can have either sign and thus 
can either act to increase or decrease the precession rate. 

The dynamical term is the one arising from the reso- 
nance and is examined in section l2~Tl The pressure term acts 
to slow the precession and it is summarised in section l2~2l 



2.1 Dynamical Precession Theory 

iHirose fc Osakil lll990l) derived the general expression for the 
ratio of the dynamical disc precession tJdyn and orbital w or b 
frequencies in terms of the mass ratio and radius of disc 
material. Their equation (8) is 
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ferumberd Il995l) is the zeroth order Laplace coefficient 
given in terms of the hypergeometric function F, q — 
M2/M1 (< 1) is the mass-ratio and r is the radius of or- 
biting material expressed as a fraction of the separation d. 
This evaluates to 
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where the coefficients are given by 
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l)Pearsonl 12003^ . iLubowl l)l992h used the fixed value of r = 
0.477 enigmatically described as "corrected for the presence 
of the companion" an d thus presumably in the limit of q — > 
0. iFrank et alJ (|l992f ). however, give the radius for j:j — 1 
resonances as, 
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This evaluates to r=0.481 for the case of j = 3 and vanishing 
q; very close to the value of the other paper but retaining 
accuracy for q 7^ 0. Substituting into © gives 
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JWarnerlll99l) . is recovered by setting j = 3 and evaluat- 
ing the summation wi th q = 0.16 . The l imiting mass ra- 
tio q ks 0.22 found by Whi tehursd l|l988a|l arises from the 
largest valu e for which rs remains with in the last stable 
stream line jMolnar fc Kobulnickvl[l992l) . Numerical simu- 
lations, however, still prod uce identifiable superhumps up 
to a mass ratio of q sa 0. 33 iWhitehurslJll994l iMurrav et alJ 
l200d) . I Whitehursd ( 1991) used an approximation for the disc 
tidal radius 



Rt ~ PRl,i 
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with (5 ~ 0.9 . When coupled to Eggleton's formula 
jEggleton II 19831) for the primary's Roche lobe radius 
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and equated to the 3:2 resonance radius in equation @, this 
gives a limiting mass ratio of g ma x = 0.28, although this is 
sensitive to the choice of f3. It should be noted, however, that 
this differs from the often cited value of g m ax = 0.33 quoted 
in that paper, as equation (5) there contains an incorrect 
power of q. A slightly less ad hoc expression fo r the tidal 
radius comes from fitting to the simulations of iPaczvriskil 
lll977l) 
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Equating this to equation @ gives a limiting mass ratio of 

q= (0.6)4j-l (15) 
which sets a maximum g ma x = 0.39 for a 3:2 resonance. 

2.2 Pressure Contribution 

ILubowl l)l992t) showed that the pressure term can be ex- 
pressed as 
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where ui p is the angular orbital frequency of a parcel of gas 
in the disc, k is the radial wavenumber of the mode and c 
is the gas sound speed. Clearly the pressure term acts in 
the opposite, retrograde sense to the dynamical term. For a 
spiral wave, the pitch angle i is related to k by 

1_ 

kr 

From the resonance condition (eg. IWarnedll995l eqn. 3.37), 
we have 
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Hence, for t he 3:2 resonance uj v = 3w or b. For the fixed radius 
r — 0.477, iMontgomervl (1200111 corrected earlier errors to 
derive a contribution 

2 
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For our general case using 11911 and where 
equation JSJ, we have 
i 

j 3 2 
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r is given by 
(21) 



vtJorbrf tani, 

To proceed further, we need to understand the behaviour of 
the final dimensionless term in brackets. Since values for c 
and i may vary according to the peculiar characteristics of 
any particular system, we will examine this in more detail 
in section l3~2l 



3 COMPARISON 

Observers normally present their measurements of the pre- 
cesion period in terms of the period excess 

Psh Porb 
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or equivalently (noting u s h = i^orb — w) 
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Since there are relatively few systems with accurately 
measured values of q, it is often convenient to make use of 
the theoretical relation 

M 2 « 0.11P O rb (26) 

llFrank et al.lll992T) that follows from the assumption that 
the secondary has a main sequence structure or an observa- 
tionally derived equivalent 

M 2 = (0.038 ± 0.003)P o 1 r b 8±0 ' 09 (27) 

JSmith fc Dhillonlll998l) . In both cases M is measured in 
solar masses and P or b in hours. 



3.1 Dynamical precession only 

For different assumed values of Mi we ca n plot theoret - 
ically predicted lines on the e-P or b plane. iMurravl j[2000jl 
used e quationElwith the fixed radius value given bv lLubowl 
( 1992) to compare the observed distribution with theory in 
just this way. He concluded that "superhumps observations 
cannot be adequately explained in terms of purely dynamical 
precession". However, including the q dependence of r given 
in equation pleads us to a different conclu sion. Figure Q re- 
produces the comparison of lMurravl 11200011 using both meth- 
ods. The lines are plotted for Mi = 0.76, 0.76 ± 0.22 and 
Mi = 1.44. This shows that the distribution is compatible 
with the boundary imposed by the condition M w( j < 1.44 
when the full q dependence is included. In fact, the most we 




Figu re 1 . Com parison of the obser ved superhump data from 
IPattersorJ ll99ft\ with the model from lMurravl fcOOttk where the 
resonant radius is independent of q (solid), and with a model in- 
cluding the q dependance (dashed). Both family of curves show 
the theoretical lines derived using the mean (0.76) and ±lcr 
(±0.22) values for Mi of all CVs and in t he M2-P or b relation 
(equation 1271 from lSmith &: DhillorJ ll998l) as well as a further 
line with Mi = 1.44. e decreases with increasing Mi. The turn 
to lower e at P or b < 1.4 may arise from the secondary becom- 
ing degenerate and deviating from a main sequence structure as 
assumed for the theoretical curves. 



can conclude is that the distribution of superhumping sys- 
tems suggests that they have a primary mass higher than the 
general CV population. This would not be entirely supris- 
ing since to be a superhumping system requires a small mass 
ratio and will thus tend to select higher Mi systems. Also 
apparent is the deviation of the secondary from a main se- 
quence structure at the short period turnoff. 

More convincing evidence for the need for a pressure re- 
lated effect on the precession rate is provided by the data for 
e and q plotte d in Figured Th e 11 sy stems on this plot are 
those used in iPatterson et alJ i2005l) . These have directly 
determined values of q. The exception is OY Car which ap- 
pears in Fig. 9 of that paper but not in the cor respondin; 
Table 7. We take the values for this system from iPattersoi 
(2001). We can see that the data are certainly not compati- 
ble with the asumption of a 3:2 resonance. If anything, they 
cluster around the prediction for a 4:3 resonance, although 
given the way the resonances 'pile up' one can argue that it 
is inevitable that some resonance would fall near the data. 
It is intriguing to note however, just how close the two most 
accurately measure values of q for the systems XZ Eri and 
DV UMa dFeline et alJl2004h lie to the 4:3 theoretical line 
(see Figure 0. It would thus be extremely interesting to 
see modelling using the method employed by these authors 
applied to the other eclipsing systems to determine simi- 
larly accurate values for q. We summarise the x 2 value for 
the data against each resonance in Table |3~T1 and note the 
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Mass Ratio (q) 

Figu re 2. Compariso n of the observed superhump data from 
iPatterson et all i2005Tl (circles) , for systems with accurately mea- 
sured values of q, and the dynamical precession rates of discs with 
radii calculated under different assumptions. The resonance pe- 
riod excesses are calculated using equation and the discs with 
radii equal to the tidal and (unphysically) the Roche lobe radius 
are calculated using equation 0. The stars mark the additional 
'challenging' systems U Gem and TV Col. 
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Mass Ratio (q) 

Figure 3. A closer view of the dat a used in the comparis on of 
the observed superhump data from IPatterson et all feOOlJl , for 
systems with accurately measured values of q, and the dynamical 
precession rates of discs at resonant radii (solid line) . Also plotted 
are the model lines including an approximation to pressure effect 
(Model A; dashed line). 







Resonance 


x 2 


3:2 


49.5 


4:3 


1.28 


5:4 


18.9 



Table 1. Summary of the x 2 statistic comparing purely dynam- 
ical resonant precession to the observed period excess. 

remarkable reduction in \ 2 f° r the 4:3 model. If a purely 
dynamical 4:3 resonance were the correct model, the proba- 
bility that the value of \ 2 would exceed the measured value 
is 0.23. 

W e have followed the precendent of IPatterson et alJ 
(2005) in the systems considered above but we mention 
here briefly two further systems: U Gem and TV Col. 
On the grounds that U Gem does not show superhumps, 
IPatterson et alJ ll2005l) used its observed mass ratio to place 
an upper limit on the system B B Dor which does have m ea- 
sured period excess. However, ISmak fc Waaeenl ll2004l) re- 
ported the detection of superhumps in the 1984 outburst of 
U Gem. As a result, we dropped BB Dor from the analy- 
sis. Similarly, iRetter et al.l ll2003l) report a superhump de- 
tection near the expected 6.3-h in TV Col from an exhaus- 
tive search in archival and from fresh 2001 observations. 
Two other campaigns, how ever, failed to confirm this re- 
sult iPatterson et alJl2005l) . Neither U Gem nor TV Col fit 
comfortably within the context of the precessing disc theory 



set out above. However, both have noticeably large errors on 
the value of their period excess. In the case of U Gem this 
reflects a systematic trend as the putative superhump period 
drifted to longer values over time. The outburst was also no- 
table for its unusual length and completely out of character 
for the system. To reach the observed period excess would 
require the disc to extend beyond the tidally truncated ra- 
dius, which is not impossible for a disc of finite viscosit y 
l|Lvnden-Bell fc Pringlelll974l: iPapaloizou fc Pringld Il977l) . 
but is uncomfortably close to the value we would get if the 
disc had a radius equal to that of the primary's entire Roche 
lobe! The most plausible explanation for such a large period 
excess and the unusual character is that the transient term 
<^tran has become important and that the change in the ob- 
served period excess reflects a change in the rate of growth 
of the disc eccentricity. The uncertainties in the values for 
TV Col make it a difficult system to assess. The expected 
disc radius is again much larger than the tidal truncation 
radius and in addition to the large range of allowed q there 
is the possibility of unknown systematic errors due to the 
multiple components of the emission lines (lHellierll993l). We 
thus n eglect this system also but echo the call of jRetter et ail 
(2003) for a search for permanent superhumps in similarly 
long period systems. 

Despite the encouraging agreement with a 4:3 reso- 
nance there are several theoretical hurdles to overcome be- 
fore we could accept such an explanation. The theory out- 
lined in Section |21 is o nly a summary of the extensive liter - 
ature in this area (eg iGoldreich fc Tremainel il978L I1979D : 
iBorderies. Goldreich fc Trern^drrel ~ (^83T) '). In outline, the 
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System 


1 


"eff 
^3 


WZ Sge 





.718 


OY Car 





.770 


XZ Eri 





.844 


IY UMa 





.787 


Z Cha 





.864 


HT Cas 





.816 


DV UMa 





.831 


OU Vir 





.762 


V2051 Oph 





.707 


DW UMa 





.873 


UU Aqr 





.886 



Table 2. Values for r B ff derive d from the observed period excesses 
for the iPatterson et alj 120051) calibration systems. 



analytical method employed is to decompose the effective 
potential into its harmonic components, to introduce this 
into the fluid equations in a standard linear perturbation 
analysis and to look at the response of the disc material. 
The behaviour of the disc in superhumping systems has a 
close analogue in galactic dynamics. The "dynamical reso- 
nance" above corresponds to the inner Lindblad resonance 
of the syste m and there is a similar res ult to that given, for 
example bv lBinnev fc Tremaind il987T) . that leads to spiral 
waves being generated in the disc. It is these spiral waves 
that couple with the tides to excite disc eccentricity. There 
is also a corotational resonance that acts to suppress eccen- 
tricity. 

Cursory consideration shows the fundamental contra- 
diction of treating the dissipation arising from an eccentric 
disc (an inherently collective phenomenon) by the preces- 
sion of single particle orbits at a resonant radius as carried 
out above. As a first attempt to correct for the untreated 
collective effects, we might assume that the precession can 
be characterised by an effective radius (r e ff) interior to r% 
that would produce dynamical precession with the observed 
period. The values for r c g d erived from the observ ed super- 
hump period excesses of the IPatterson e t al. (2005) systems 
with well measured q are given in Table loTTl Assuming that 
the ratio r e g/r3 is a constant for all systems, we derive a 
best value for it of 0.827. This places r e g extremely close to 
r4 with a barely different value for the total \ 2 - in fact, the 
radius calculated from a j = 4 resonance differs from this 
best possible radius by a suprisingly small 0.2%! 

The above result not withstanding, a treatment that 
deals explictly with the coupling of the tides, spiral arms 
and disc eccentricity explicitly ought to be p referred. This 
is exactly the approach used bv lLubowl Jl991al) to derive the 
additional term related to pressure effects which we turn to 
next. One of the important results from that paper was the 
recognition that the 3:2 resonance "is unique in that it is the 
innermost resonance for which an eccentric Lindblad reso- 
nance appears without an overlapping eccentric corotational 
resonance. . . This property allows that resonance to easily 
excite eccentricity." Such a conclusion is a strong argument 
against any dynamical resonance other than 3:2 being an 
important factor in the disc behaviour. The ability of a sec- 
ondary magnetic field to give rise to a precession rate charac- 
teristic of a hi gher resonance dPearson. Wvnn fc Kine|ll997t 
|Pearsoni r2003) may reflect the fact that the perturbation in 



that case is no longer small and capable of represention by 
a linear analysis. 



3.2 Inclusion of Pressure Effects 

Proceeding under the assumption that the excited resonance 
is 3:2 and that the difference between the measured and 
expected e is due to the pressure effect we can update the 
analysis of lMurravl |2000) for all the systems with measured 
q. These are summarised in tablc !3.2l Fo r the final col umn we 
have used values for Mi from IPat terson et alj <|200 5T) exce pt 
for V2051 Oph which we take from lRitter fc Kolbl hoO^i . 

To make progress with our comparison we need val- 
ues for c and i in equation 1211 We will consider two cases: 
A) where ° co ^ is fixed for all systems and B) where c is 
evaluated using analytic expression taken from detailed disc 
models. 

3.2.1 Model A 

If we assume that the final bracketed term of equation 12H 
is constant we can rewrite it as 

^=-^(1 + ^ (28) 
where 

»-\&- < 29 » 

iLubowl l(l"992T) gave a range of 0.01-0.05 for the ratio — —, -,. 
iMurravl (|2000) used a value of 0.05 in his considerations to 
derive values for i. Using th is value and the me an i = 17° 
found from simulations by iMontgomervl l|2QQlf), we have 
7/a = 0.0134. Compared to the data for the eclipsing sys- 

2 

terns, this gives — = 3.3. Allowing tia to be a free param- 
eter, however, we derive an optimal value of 7/a = 0.0107 
with a corresponding N * _ 1 = 1.42. Figure |31 shows the 
comparison for this value in graphical form. 

The corresponding plot to Figure^ comparing a purely 
dynamical 3:2 resonance and a model including pressure, is 
shown in Figure 01 It should be noted here that the val- 
ues assumed for Mi in p roducing the lines come from the 
weighted mean derived bv lSmith fc Dhillonl Jl998l) for dwarf 
nova systems ( Mi = 0.69 ± 0. 01) rather than that for all 
systems used bv IMurravl <|2000F) . W e have also updated th e 
data to the latest compilation from IPatterson et alJ i2005^ . 
Coupling the mean primary mass with q max = 0.39, we can 
derive an equivalent P or b,max = 3.5-h. For the ultimate re- 
striction Mi < 1.44M© we have P rb,max = 5.5-h (coinci- 
dentally the same as that of TV Col). 

We might like to compare t he implied distribution o f 
Mi from our model with that of ISmith fc Dhillonl Jl998i) . 
However, this is precluded by the selection effect alluded to 
in section 13.11 being at work. Systems with longer periods 
(and thus higher M2) can accomodate more massive pri- 
maries and still fit within the g max limitation. Hence, the 
distribution of Aii would be expected to differ between all 
dwarf novae and the superhumping subset. Figure 0] a ls° 
shows systems lieing above the limiting e allowed by our 
formulation of the pressure term. This, along with the range 
of derived t/a shown in Table IqT21 and the poorer fit to the 
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System 




<^dyn 


^prcss 




c 

ui alb d 


c 




(d- 1 ) 


(d" 1 ) 


(d- 1 ) 






(10 4 m s" 1 ) 


WZ Sge 


1.01 


2.13 


-1.12 


0.12 


0.036 


2.07 


OY Car 


1.98 


3.59 


-1.61 


0.14 


0.044 


2.16 


XZ Eri 


2.70 


4.02 


-1.32 


0.13 


0.039 


2.02 


IY UMa 


2.15 


3.71 


-1.56 


0.15 


0.047 


2.29 


Z Cha 


2.96 


4.17 


-1.21 


0.13 


0.041 


1.80 


HT Cas 


2.73 


4.34 


-1.62 


0.15 


0.047 


2.14 


DV UMa 


2.43 


3.73 


-1.30 


0.15 


0.046 


2.38 


OU Vir 


2.73 


4.99 


-2.26 


0.18 


0.055 


2.87 


V2051 Oph 


2.93 


6.20 


-3.27 


0.20 


0.061 


3.21 


DW UMa 


2.78 


3.79 


-1.01 


0.16 


0.049 


1.99 


UU Aqr 


2.52 


3.33 


-0.81 


0.16 


0.048 


1.78 



Table 3. Derived pressure force contribution to the precession of the lPatterson et al.1 120051) calibration systems. The last two columns 
assume i = 17° . 




orb 



(hrs) 



Figu re 4. Comparis on of the observed superhump data from 
iPatterson et all 120051) with a model with purely dynamical pre- 
cession (solid) and including the pressure effect (dashed). Both 
family of curves show the theoretical lines derived using the mean 
(0.69) and ±1ct (0.01) values for Mi foun d for dwarf novae and in 
the M 2 -P olb relation (equation^ from lSmith fc Dhilloij Il998f) 
as well as a further line with Mi = 1.44. e decreases with increas- 
ing Mi. The horizontal lines mark the limiting precession rate 
derived from g max = 0.39 in either case. 



data than a simple resonance model, suggests that a one size 
fits all value for r/ is not appropriate. 



3.2.2 Model B 

Since the effect of the pressure term relies on the sound speed 
in the disc, we turn to detailed models of hot discs to eval- 
uate this i n terms of fundamental parameters. From equa- 
tion Al of ICannizzo fc Reifll lll992T) we have the mid-plane 
temperature in terms of the mass transport rate through the 
disc 



m /64 a \-ttj /fi H \4 i i 



mV 

2^ ) 



(30) 



where a is the Stefan-Boltzman constant, R is the gas con- 
stant, qh and /ih are the Sunyaev-Shakura viscosity param- 
eter and mean molecular weight in the hot state respectively 
and assuming an opacity k = KppT~ 3 5 where kq = 2.8 x 
10 23 m 2 kg -1 is appropriate (ICannizzo. Shaffer fc Wheeler I 
1988). We can find a suitable value o f M for an outb urst- 
ing dwarf nova using the approach o f ICannizzo! <ll993l) and 
ICannizzo. Shafter fc Wheeler I Jl988l) by assuming that the 
disc fills to a mass /Af max . Af max is the maximum mass the 
disc can hold in the cold state without exceeding E max at 
some point ie. 



fM n 



= / 



2-7rr£ max (r)dr 



f 27rr| 
J 3.05 n 



c(rd). 



Now, equation A3 of ICannizzo fc Reifll il99: 

face density in the hot state as 



E = 



/64_o_\ io / /hh 
V 9 no) V R 



405 kg m /i*a H B M- 



M 



n 



10 8 m 



10- 



M y 



(31) 
(32) 

gives the sur- 
(33) 



(34) 



where we have used u p = (GMiMqt -3 ^" 1 . Integrating this 
and equating it to the expression for /M max from 1)32(1 . we 
can rearrange for the mass transport rate through the disc 



M = 9.67 x 10" 



kgs IotJ Km) ^ 



xr 2 - 57 ' 



which, with equation (IMOfl . gives 
2 jkT 



(35) 



(36) 



/ Qh 

Vo.i 



(— V 

1 0.02/ 
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w -1.071 0.771 ,,-0.258 



\oaJ 



(37) 



Finally, we can combine this with equations @, pfip . I|17|l 
and 1191 1 to get 

-2.83 x 10 + cot i f ^j-J 



I^orb 



QC 
02 

x 

0.4 



M -1.071 r 0.771 M -0.258 



(38) 



Since we want an expression for impress in terms of q only, we 
look to eliminate w or b b y using a form of the m ass-radius 
relation recommended bv lSmith fc Dhillonl il998f) 



^ = (0.91 ± 0.09)M 2 075±0 ' 04 . (39) 

When this is equated to the size of the secondary's Roche 
lobe from 113L we can rearrange for the separation 



d 



O.91M 2 °- 75 i? 



E(q) ' 
Using this with Kepler's Law 

to arh = [GMi(l + q)M e ]i d~i 



(40) 



(41) 




and a disc radius = /3J?l,i, we arrive at the final expres- 



Mass Ratio (q) 

Figu re 5. C omparison of the calibration systems from 
iPatterson et al] 120051) with a model with purely dynamical pre- 
cession (solid) and including the detail pressure term in equa- 
tion g2J (dashed). 



^orb 



where 



I IB 



-J°TJB- 



[Eiq- 1 ) 



[E(q)Y 



(1 + 9)' 



(42) 



0.0209 



cot^r 42 ^ 

\0.lJ V0.02 







(43) 



Hence, we can see that although the pressure term is not 
purely expressible as a function of q alone, the various pa- 
rameters are either expected to remain reasonably constant 
between systems (eg. jj,n) or enter in with a weak depen- 
dency such as Mj -0 ' 242 . As mentioned above, the observed 
distribution for all CVs is Mi = 0.76 ± 0.22 which would 
produce a variation of only ~ 7% in the predicted pressure 
contribution. The dwarf nova subsample had an even smaller 
range for Mi. 

The final expression for the precession rate w as fitted 
to the calibration systems of IPatterson et all i2005t) using a 
single free parameter as the constant of proportionality to 
the functional form of q in (I42H . The best fitting value of 

2 

7/b = 0.0109 gives N x _ 1 = 1.04. With the typical values 
used above this implies i = 61°. A comparison with the data 
is plotted in Figure For the lowest q systems the effect of 
the pressure term actually makes e negative, ie. force the 
precession to become retrograde. 

Inverting the process for systems of known e but un- 
known q a llows us to derive valu es for the 88 systems in 
Table 9 of IPatterson et alJ (120051) . These are listed in Ta- 
ble !3.2.2l The errors quoted in the table reflect the observa- 
tional errors propogated as appropriate but not the system- 
atic e rrors that may ari s e from variations of the parameters 
m the lSmith fc Dhillonl lll998l) relations. It is apparent that 



some systems have predicted values of q in excess of the ex- 
pected limit of g ma x = 0.28 for (5 — 0.9. The largest derived 
value of q = 0.437 can be accommodated if /3 — 0.94. If 
the resonance were to transition to 4:3 at some g max , the 
derived values of q would be even higher. However, axisym- 
metric structure models were used in the above derivation. 
As Vj approaches _Rt this will become an increasingly invalid 
assumption. It might be reasonably expected then that these 
high q systems have true mass ratios close to g max although, 
as we have seen, this limit is not well determined. 

It is also worth noting just how close is the derived Mi 
in EG Cnc to the Chandrasekhar limit. Although the mass- 
period relation we have used ('eauation l391 was derived from 
observations, there are strong theoretical grounds for expect- 
ing a significantly different relation for the lowest mass sec- 
ondaries once they become degenerate or semi-degenerate. 
EG Cnc is towards the low-end in the range of M 2 where this 
effect may be important and we may be extrapolating the 
relation beyond its range of validity. A non-main-sequence 
star would be expected to have a larger radius for the same 
mass as a main-sequence equivalent. This would lead to a 
larger pressure effect for a given value of q. As a result we 
would tend to underestimate both M 2 and q when such de- 
viation becomes important. In principle, analytic forms for 
the mass-radius relations for the appropriate ranges of Mi 
could be included in place of 1391 . 

As a disc outburst proceeds, we would expect M to 
steadily decrease. Since c « Mu this would cause the 
pressure term to also shrink and thus the period excess to 
increase during an outburst. Observationally, the opposite 
appears t o be the case, with th e period excess decreasing 
with time (IPatterson et alll995t) . The only available param- 
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Name 


Porb(d) 


e 










M 2 




Mi 


DI UMa 


0.05456 


1) 


0.0133 


( 6) 





071 


: 2) 





058 





821(25 


V844Her 


0.05464 


1) 


0.0243 


( 9) 





111 


; 3) 





058 





526(16 


LLAND 


0.05505 


1) 


0.0290 


36) 





128 


14) 





059 





459(49 


SDSS0137-09 


0.05537 


4) 


0.0248 


20) 





113 


: 7) 





060 





528(35 


ASAS0025+12 


0.05605 


5) 


0.0206 


21) 





097 


: s) 





061 





624(49 


AL Com 


0.05667 


3) 


0.0120 


( 7) 





066 


; 3) 





062 





933(35 


WZ Sge 


0.05669 


1) 


0.0092 


( 7) 





056 


: 3) 





062 


1 


100(49 


RX1839+26 


0.05669 


5) 


0.0173 


20) 





085 


7) 





062 





724(61 


PU CMa 


0.05669 


5) 


0.0222 


20) 





103 


: 7) 





062 





599(43 


SW UMa 


0.05681( 


14) 


0.0245 


27) 





112 


10) 





062 





556(50 


HV Vir 


0.05707 


1) 


0.0200 


( 9) 





095 


; 3) 





062 





657(23 


MM Hya 


0.05759 


1) 


0.0184 


10) 





089 


14) 





063 





710(29 


WX Cet 


0.05829 


4) 


0.0199 


15) 





095 


: 5) 





065 





682(39 


KV Dra 


0.05876 


7) 


0.0233 


22) 





107 


: s) 





065 





610(46 


T Leo 


0.05882 


1) 


0.0236 


14) 





108 


: 5) 





066 





605(29 


EG Cnc 


0.05997 


9) 


0.0067 


( 8) 





047 


: 3) 





068 


1 


433(88 


V1040 Cen 


0.06028( 


10) 


0.0310 


27) 





136( 


10) 





068 





501(38 


RX Vol 


0.06030( 


20) 


0.0178 


20) 





087 


: 7) 





068 





782(65 


AQ Eri 


0.06094 


6) 


0.0284 


21) 





126 


; s) 





069 





549(34 


XZ Eri 


0.06116 


1) 


0.0270 


16) 





121 


: e) 





070 





576(28 


CP Pup 


0.06145 


6) 


0.0171 


20) 





085 


: 7) 





070 





830(71 


V1159 Ori 


0.06218 


1) 


0.0320 


11) 





140 


;4) 





072 





512(15 


V2051 Oph 


0.06243 


1) 


0.0281 


25) 





125 


9) 





072 





576(43 


V436 Cen 


0.06250( 


20) 


0.0212 


32) 





099 


12) 





072 





725(85 


BC UMa 


0.06261 


1) 


0.0306 


14) 





134 


: 5) 





072 





538(21 


HO Del 


0.06266( 


16) 


0.0276 


35) 





123 


13) 





072 





588(63 


EK TrA 


0.06288 


5) 


0.0321 


25) 





140 


10) 





073 





519(35 


TV Crv 


0.06290( 


20) 


0.0325 


32) 





142 


12) 





073 





513(44 


VY Aqr 


0.06309 


4) 


0.0203 


15) 





096 


: 5) 





073 





761(43 


OY Car 


0.06312 


1) 


0.0203 


15) 





096 


; 5) 





073 





761(43 


RX1131+43 


0.06331 


8) 


0.0259 


16) 





117 


: e) 





074 





630(32 


ER UMa 


0.06336 


3) 


0.0314 


11) 





138 


;4) 





074 





536(17 


DM Lyr 


0.06546 


6) 


0.0281 


31) 





125 


12) 





078 





620(58 


UV Per 


0.06489( 


11) 


0.0234 


23) 





108 


; s) 





077 





711(56 


AK Cnc 


0.06510( 


20) 


0.0368 


33) 





158( 


13) 





077 





486(40 


AO Oct 


0.06557( 


13) 


0.0242 


39) 





111 


14) 





078 





704(92 


SX LMi 


0.06717( 


11) 


0.0347 


25) 





150 


10) 





081 





538(35 


SS UMi 


0.06778 


4) 


0.0360 


15) 





155 


: e) 





082 





528(20 


KS UMa 


0.06796( 


10) 


0.0241 


30) 





110 


ii) 





082 





747(75 


V1208 Tau 


0.06810( 


20) 


0.0374 


28) 





161( 


ii) 





083 





514(35 


RZ Sge 


0.06828 


2) 


0.0306 


28) 





134( 


ii) 





083 





617(49 


TY Psc 


0.06833 


5) 


0.0347 


15) 





150 


; e) 





083 





553(21 


IR Gem 


0.06840( 


30) 


0.0351 


66) 





152( 


26) 





083 





548(93 


V699 Oph 


0.06890( 


20) 


0.0197 


28) 





094( 


10) 





084 





895(97 


CY UMa 


0.06957 


4) 


0.0364 


14) 





157 


: 5) 





085 





545(19 


FO And 


0.07161( 


18) 


0.0349 


40) 





151 


16) 





089 





592(61 


OU Vir 


0.07271 


1) 


0.0326 


15) 





142 


: e) 





092 





643(26 


VZ Pyx 


0.07332 


3) 


0.0333 


20) 





145 


;s) 





093 





641(34 


CC Cnc 


0.07352 


5) 


0.0487 


27) 





207 


ii) 





093 





449(25 


HT Cas 


0.07365 


1) 


0.0330 


30) 





144( 


12) 





093 





651(52 


IY UMa 


0.07391 


1) 


0.0260 


13) 





117 


;s) 





094 





802(33 


VW Hyi 


0.07427 


1) 


0.0331 


(8) 





144 


;s) 





095 





658(13 


Z Cha 


0.07450 


1) 


0.0364 


(9) 





157 


;4) 





095 





607(14 


QW Ser 


0.07453( 


10) 


0.0331 


40) 





144( 


15) 





095 





661(71 


WX Hyi 


0.07481 


1) 


0.0346 


14) 





150 







096 





639(23 


BK Lyn 


0.07498 


5) 


0.0479 


(7) 





204 


;s) 





096 





471( 7 


RZ Leo 


0.07604 


1) 


0.0347 


25) 





150( 


10) 





098 





654(42 


AW Gem 


0.07621( 


10) 


0.0422 


27) 





180 


11) 





099 





548(33 


SU UMa 


0.07635 


5) 


0.0317 


12) 





139 







099 





713(23 



The errors on q reflect the 
No allowance has been made 



Table 4. Derived parameters for all the systems with measured period excesses from lPatterson et alj 1200. 
observational errors propogated appropriately. M2 has been calculated using equation 1271 and Mi = 
for the systematic error that would arise from the errors in the parameters in this relation. Similarly, no error is quoted for M.% since this 
is dominated by the assumption of main sequence structure. 
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Star 


Porb(d) 


e 




M 2 


Mi 


SDSS1730+62 


0.07655( 9) 


0.0376(22) 


0.162( 9) 





099 


0.615(33) 


HS Vir 


0.07690(20) 


0.0477(23) 


0.203(10) 





100 


0.493(23) 


V503 Cyg 


0.07770(20) 


0.0430(27) 


0.183(11) 





102 


0.555(34) 


V359 Cen 


0.07990(30) 


0.0388(40) 


0.166(16) 





106 


0.639(61) 


CU Vel 


0.07850(20) 


0.0293(36) 


0.130(14) 





103 


0.798(83) 


NSV 9923 


0.07910(20) 


0.0412(30) 


0.176(12) 





105 


0.595(41) 


BR Lup 


0.07950(20) 


0.0340(40) 


0.147(15) 





105 


0.716(75) 


V1974 Cyg 


0.08126( 1) 


0.0471(10) 


0.201( 4) 





109 


0.544(11) 


TU Crt 


0.08209( 9) 


0.0397(22) 


0.170( 9) 





111 


0.653(34) 


TY PsA 


0.08414(18) 


0.0417(22) 


0.178( 9) 





115 


0.648(32) 


KK Tel 


0.08453(21) 


0.0368(31) 


0.158(12) 





116 


0.734(56) 


V452 Cas 


0.08460(20) 


0.0497(33) 


0.212(14) 





116 


0.550(37) 


DV Uma 


0.08585( 1) 


0.0343(11) 


0.149( 4) 





119 


0.801(23) 


YZ Cnc 


0.08680(20) 


0.0553(26) 


0.236(12) 





121 


0.514(25) 


GX Cas 


0.08902(16) 


0.0449(25) 


0.191(11) 





126 


0.660(37) 


NY Ser 


0.09775(19) 


0.0623(35) 


0.268(16) 





146 


0.546(33) 


V348 Pup 


0.10184( 1) 


0.0640(40) 


0.276(19) 





156 


0.565(39) 


V795 Her 


0.10826( 1) 


0.0760(10) 


0.336( 5) 





172 


0.512( 8) 


V592 Cas 


0.11506( 1) 


0.0625( 5) 


0.269( 2) 





189 


0.703( 6) 


TU Men 


0.11720(20) 


0.0717(32) 


0.314(16) 





195 


0.621(32) 


AH Men 


0.12721( 6) 


0.0887(16) 


0.406( 9) 





222 


0.546(13) 


DW Uma 


0.13661( 1) 


0.0644(20) 


0.278(10) 





248 


0.893(31) 


TT Ari 


0.13755( 1) 


0.0847( 7) 


0.383( 4) 





251 


0.655( 7) 


V603 Aql 


0.13810(20) 


0.0572(51) 


0.244(23) 





252 


1.032(97) 


PX And 


0.14635( 1) 


0.0898(14) 


0.412( 8) 





277 


0.671(13) 


V533 Her 


0.14730(20) 


0.0719(20) 


0.315(10) 





279 


0.888(29) 


BB Dor 


0.14920(10) 


0.0939(10) 


0.437( 6) 





285 


0.652( 9) 


BH Lyn 


0.15575( 1) 


0.0790(30) 


0.352(16) 





305 


0.868(40) 


UU Aqr 


0.16358( 1) 


0.0702(14) 


0.306( 7) 





330 


1.077(25) 



Table 4. Cont. Derived parameters for all the systems with measured period excesses from lPatterson et alj 1200. 




P orb ( h ^) 

Figu re 6. Compariso n of the observed superhump data from 
IPatterson et al.l J2005f) with a model including the full pressure 
effect. The curves use the mean (0.69) and ±1<t (0.01) values 
for Mi found for dwarf novae and in the M2-P rb relation from 
ISmith fc Dhillonl ( 1998) as well as a further line with Mi = 1.44. 
e decreases with increasing Mi. 



eter to counter this is tani implying that the pitch angle i 

increases during an outburst. 

Although the treatment of lLubowl lll991al) includes the 
coupling of the tides, spiral waves and eccentricity, it shares 
the shortcoming of a purely dynamical resonance in that the 
properties are characterised by those at a single radius. A re- 
cent preprint of a paper bv lGoodchild fc Qgilviei <l2006l) has 
addressed this problem with a detailed analysis solving the 
equations for the disc behaviour integrated over the whole 
range of disc radii. Their final equation describing the gen- 
eration, damping and dynamics of eccentricity in the disc 
is 

=flf + ^^((7-- b )^||) + 
^ (&| (§) E) + 2^rQE6{r - r rcs ). (44) 

Here E is the complex eccentricity E = ee tu} , p and p are the 
local disc pressure and density, 7— icvb is a complex adiabatic 
exponent, £ = 2.08aj or bq ,2 r r cs is the eccentricity gr owth rate 
for a resonance at radius r rcs from iLubowl Ill991ar) and VL is 
the angular velocity of the orbiting disc material. Unsupris- 
ingly, the full solution has to be found numerically. These 
authors do so with undisturbed (vertically integrated) struc- 
ture distributions given by 
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where r ou t and r- m are the outer and inner radii disc radii and 
P sc and E sc are the pressure and surface density at the scal- 
ing radi us r sc . Considering the fir st bracketed terms in either 
case, the lCannizzo fc ReifJ (I1992T) equations used above have 
the same radial dependence. The second bracketed terms re- 
flect the behaviour close to the inner radius. The final terms 
were chosen to implement the boundary conditions at the 
outer edge of the disc. Since we would expect the inner and 
outer terms to only become important close to the limiting 
mass ratios and given the similarity otherwise of the radial 
dependences of the equations to those used earlier, we do 
not need to repeat the integration here but compare the 
final results produced. 

The results in Figure 8 of lGoodchild fc Qgilviel (l2006r) 
show a similar excursion to a negative period excess for small 
q. The curves also turn up to high e as they assymptotically 
approach gmax which our expression does not produce. How- 
ever, the turn off occurs very close to g ma x for the best fitting 
curve and the details of the behaviour here depend on the 
functional form chosen to implement the boundary condi- 
tions for the outer edge of the disc in equations H45H and 
14611 . Away from the extreme mass ratios the two sets of 
results are very close. 



3.3 Empirical Fit 

In an attempt to calib r ate a n empirical relation between e 
and q, iPatterson et alJ (120051) extended the earlier fit of 



e = 0.22g 

JPattersonll200lfl to 
e = 0.18g + 0.29g 2 . 



(47) 



(48) 



We can view these as phenomenologically derived equiva- 
lents to the Maclaurin series for our analytic expression. 
Given the complexity of the full expressions, however, rather 
than carry out the necessary differentiation, it is simplest to 
generate e(q) predictions from our formula numerically and 
then find the best fitting polynomial to these values using 
standard methods iPress et al.lll992T) . This approach gives 
an approximate formula 



e = 3.5 x 10" 4 + 0.24g - 0.12g 



(49) 



over the range 0.01 < q < 0.4. This expression differs both 
with a very small constant offset and with a quadratic term 
that causes a curvature in the opposite sense to the empirical 
expression. This can be attributed to r emoving the q ma x lim- 
itation imposed by U Gem for BB Dor. lGoodchild fc Qgilviel 
(2006) derive a best fit formula for their numerical integra- 
tions of lEBl of 



-4.1 x 10~ 4 + 0.2076q. 



(50) 



These polynomial forms are compared graphically in Fig- 
ure |7| The linear empirical fit, the full integration polyno- 
mial and our result all appear in good agreement. 



4 SUMMARY 

We have shown that the standard dynamical method of cal- 
culating the precession rate of superhumping CVs with a 




Mass Ratio (q) 

Figur e 7. Compariso n of the observed superhump data from 
IPatterson et al.l i2005t) with the polynomial approximations for 
empirical fits (equations 1471 and 1481 (solid), this work 1491 
(dashed) and full integration I5UI (dot-dashed). 



3:2 provides such a poor fit to the data that a 4:3 reso- 
nance is actually a better fit! We have confirmed the im- 
portance of including the pres sure related te rm in the cal- 
culation of the precession rate jMurravll2000F) which fits the 
data better than any pure resonance model. The pressure 
term has been reduced to a function of q and the total, an- 
alytic precession rate shown to be equivalent to the empir- 
ically derived expressions of IPatterson et all i2005tl . These 
anayltic expressions also produce precession rates in good 
agreem ent with those from the deta iled integrations carried 
out by iGoodchild fc Ogilviel teOOrJ) . This formulation has 
been used to calculate values for q in systems which would 
otherwise be unknown. 
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